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The spin dynamics of bilayer cuprate compounds are studied in a basic model. The magnetic 
spectral properties are calculated in linear spin-wave theory for several stripe configurations which 
differ by the relative location of the stripes in the layers. We focus on the bilayer splitting of 
the magnon bands near the incommensurate low energy peaks as well as near the ir resonance, 
distinguishing between the odd and even channel. We find that a x-shaped dispersion near the 
7r resonance is generic for stripes. By comparison of our results to neutron scattering data for 
YBa2Cu306+ x we conclude that the stripe model is consistent with characteristic features of bilayer 
high-Tc compounds. 
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I. INTRODUCTION 

Subsequent to predictions of stripe formation; 1 : 2 : 3 
characteristic signatures of spin- and charge order 
have been found in a variety of high-T c cuprate 
superconductors, including La 2 - x Sr x Cu04 (LSCO) 
and YBa 2 Cu306+ x (YBCO). Neutron scattering 
experiments^ have provided evidence for spin order at 
low energies through a pattern of incommensurate peaks 
around the antiferromagnetic wave vector. Although 
more difficult to detect, charge order has been observed 
in LSCO co-doped with Nd 6 - as well as in YBCO without 
codopingi 

Since LSCO and YBCO are paradigmatic for mono- 
layer and bilayer compounds, stripe-like "low" -energy 
response is characteristic for both classes of materi- 
als. On the other hand, at "high" energies spin fluc- 
tuations appeared to be qualitatively different since a 
commensurate 7r-resonance had been observed only in 
bilayer compounds, notably in YBa 2 Cu306+x^' and 
Bi 2 Sr 2 CaCu208+xr^ whereas it seemed to be absent in 
monolayer compounds. This apparent distinction be- 
tween mono- and bilayer compounds lost its justification 
only recently, when the 7r-resonance was discovered in 
Tl 2 Ba 2 Cu06+ x , ii as the first monolayer compound. The 
fact that the 7r-resonant mode has not been detected in 
LSCO so far can possibly be ascribed to a larger effective 
strength of disorder, since the Sr-dopants are randomly 
distributed whereas in the oxygen doped compounds the 
access oxygen orders in chains. Thus, one may believe 
that, in principle, mono- and bilayer compounds have 
qualitatively similar features also at higher energies. This 
universality of low- and high-energy features calls for an 
even more unifying framework^ 

In a recent article^ we have analyzed an elementary 
monolayer model assuming that charges form a perfectly 
ordered site-centered stripe array which imposes a static 
spatial modulation of spin-exchange couplings. The re- 
sulting spin dynamics was studied using linear spin- wave 
theory. As a result, we found that the incommensu- 
rability and the 7r-resonance appear as complementary 
features of the band structure at different energy scales. 



Furthermore, the doping dependence of the resonance fre- 
quency was found in good agreement with experimental 
observations. 

In this work we extend this model to bilayer systems 
in order to predict the corresponding features of the 
magnon band structure and the magnetic structure fac- 
tor. Within each layer, holes are assumed to form unidi- 
rectional site-centered stripes. We consider several possi- 
bilities (parallel and perpendicular relative orientations) 
of the charge order in the antiferromagnetically coupled 
neighboring layers. The band structure and the T = 
inelastic structure factor for even and odd excitations are 
calculated in linear spin-wave theory. Particular atten- 
tion is paid to the band-splitting in the vicinity of the 
antiferromagnetic wave vector and to the influence of the 
interlayer coupling on the 7r-resonance energy. 

The outline of this paper is as follows. In Sec. [Q] 
the spin-only model for a bilayer system is introduced 
and motivated. Classical ground states and the result- 
ing phase diagrams for competing types of magnetic or- 
der are obtained. They are needed as starting point for 
the linear spin-wave theory. A customized formulation 
thereof is outlined in Sec. IIIII The results, namely the 
spin- wave band structure, the zero-temperature structure 
factor for even and odd excitations, and the dependence 
of the band splitting at the antiferromagnetic wave- vector 
on the strength of the interlayer coupling are presented 
in Sec. IIVI and compared to experiments in Sec. IV1 



II. MODEL 

Stripes are a combined charge- and spin-density wave. 
If the charge period pa is a multiple of the Cu spacing a 
with integer p, lock-in effects tend to suppress phason-like 
fluctuations of the density modulation. In a reductionist 
real-space picture, one may think of the holes forming 
parallel site-centered rivers of width a, which act as an- 
tiphase domain boundaries for the antiferromagnetic spin 
domains in between. 6 This implies that the period of the 
spin modulation is twice that of the charge modulation. 

To implement that the charge stripes act like antiphase 
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boundaries we follow our previous worki£ and choose the 
simplest possible implementation of exchange couplings 
within the layers stabilizing this magnetic structure: an- 
tiferromagnetic exchange couplings J between neighbor- 
ing spins within the domains and antiferromagnetic cou- 
plings A J between closest spins across a stripe. 

In our previous workli we have studied this model for 
a single layer allowing for diagonal and vertical stripe ori- 
entations. Here we focus on vertical stripes as observed 
in the superconducting cuprates and restrict our analysis 
to the representative case p = 4. This corresponds to a 
doping of one hole per 8 Cu sites since the rivers have 
a line charge of only half a hole per lattice constant. In 
addition to the in-plane couplings we consider an antifer- 
romagnetic exchange pj between two layers (cf. Fig.^l. 
The Hamiltonian of this bilayer model is given by 



H 



Ha + #1,2) 



a=l,2 



H a = i^J Q (r,r')S a (r)S a (r'), 

r,r' 

H h2 = ^^Si(r)S 2 (r), 



(la) 
(lb) 
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where r specifies the square-lattice position and a = 1, 2 
numbers the layers. The asterisks indicate that the sums 
do not include positions of charge rivers. The in-plane 
couplings J Q (r, r') defined in the text above are illus- 
trated in Fig. ^ They explicitly depend on the layer 
index if the charge distribution is different in both lay- 
ers. 

For simplicity we neglect spin anisotropy, the weak 
3D coupling between bilayers, and more compli- 
cated exchange processes such as cyclic exchange or 
Dzyaloshinskii-Moriya interactions, which all may be im- 
portant for quantitative purposes. Obviously, this simple 
spin-only model does not account for electronic correla- 
tion effects, e.g. a spin gap at low energies due to the 
formation of Cooper pairs is not incorporated. Never- 
theless we expect that our model provides a qualitatively 
adequate description of the spin fluctuations well above 
the gap energy. 

The actual stripe configuration is determined by sev- 
eral influences. Besides the magnetic exchange energy 
one also has to take into account the Coulomb energy, 
and in principle also a further reduction of the four- 
fold symmetry of Cu02 planes in orthorhombic struc- 
tures which may favor a certain alignment of the stripes. 
In YBa2Cu306+ x the formation of CuO chains along the 
6-direction may favor a parallel alignment of stripes. 

We find that three different stripe configurations may 
be realized physically (see Fig. ^) . The exchange energy 
favors parallel stripes lying exactly on top of each other. 
This configuration is free of magnetic exchange frustra- 
tion, each bond can be fully saturated. However, this 
configuration is disfavored by the Coulomb energy which 
would favor a configuration where stripes are parallel but 




FIG. 1: Classical ground states for bilayer systems with par- 
allel (upper row), shifted parallel (middle row), and perpen- 
dicular hole stripes (lower row) for a stripe spacing p = 4. 
The exchange couplings of the simple model are illustrated 
in the lower row: AF couplings J > for nearest neigh- 
bors within the domains (bold dashed), A J between nearest 
neighbors across a hole stripe (zig-zag) and couplings pj be- 
tween spins one above the other (dashed). Frustration of ex- 
change coupling may lead to a canting of spins (calculated for 
p — 0.09 and A = 0.07 in the middle and bottom row, respec- 
tively). Possible magnetic unit cells are outlined by gray lines, 
identical gray levels of spins correspond to identical canting 
angles. 



shifted with respect to each other by half a stripe spac- 
ing. (In our simple model, where holes are assumed to be 
site centered, this configuration is only compatible with 
even stripe spacings p.) The gain in Coulomb energy has 
to be paid by a loss of exchange energy. For certain pa- 
rameters, a third configuration may be favorable, where 
the charge stripes of the two layers are perpendicular. 

For the later analysis it is instructive to anticipate 
that for these configurations the Hamiltonian has discrete 
symmetries. We focus on symmetries involving an ex- 
change of layers. For parallel and shifted parallel stripes, 
this symmetry is just the reflection z — > —z combined 
with a translation (coordinates are chosen such that the 
planes are parallel to the xy plane). For perpendicular 
stripes, one needs to add a rotation around the z axis. 



A. Energetic estimates 

To estimate the Coulomb energy for the three stripe 
configurations, we assume a charge-density modula- 
tion p(r) = pi(r)8(z) + p 2 {r)5(z - d) with p a (r) = 
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pa^ cos(k a r) where the planes separated by d are per- 
pendicular to the z-direction. For simplicity, only the 
first harmonic of the charge modulation is retained. Par- 
allel stripes are described by ki = k2 = ke x and 
pf^ = p^ = p, shifted parallel stripes are realized for 
ki = k2 = ke x and p^' = —p% — P, and perpendicular 
stripes for ki = ke x , k2 = ke y and pf 1 — p% = p. For a 
stripe spacing pa the charge- modulation wave vectors are 
given by k = 2ir/(pa), the amplitude by p — e/(2pa 2 ). 
Calculating the Coulomb coupling energy per square lat- 
tice site 
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where A denotes the area of the planes, we find, in the 
limit A — > oo, a vanishing Coulomb coupling for perpen- 
dicular stripes, an energy cost 
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for parallel stripes, and an energy gain of the same size 
for shifted parallel stripes. For YBCO with a « 3.85 A, 
d « 3.34 A, J = 125 meV, S = i and for a stripe spacing 
p = 4 we obtain AEc ~ 29 meV. 

For antiferromagnetic YBCO the magnetic interlayer 
superexchange is reported to be p ss 0.08»i* For paral- 
lel stripes, spins are not frustrated and, in a classical 
picture, antiparallel in different layers, Si(r) = — S2(r). 
Thus, the exchange coupling roughly leads to an energy 
gain of order pJS 2 ~ 3 meV, whereas the energy gain 
will be smaller for the other two configurations due to 
frustration. 

Thus, within our rough estimate, the Coulomb energy 
appears to be up to one order of magnitude larger than 
the exchange energy, such that one might expect the par- 
allel shifted configuration to be the only physical one. On 
the other hand the actual Coulomb energy may be sig- 
nificantly smaller than the result of our estimate since 
we have completely neglected screening. For almost un- 
doped YBCO a relative large value of e m 15 for the 
static dielectric constant at T = 4K is reported*^ There- 
fore the Coulomb energy might be of the same order of 
magnitude as the magnetic exchange energy. Due to the 
crudeness of our estimate no stripe configuration can be 
strictly ruled out. 



B. Classical ground states 

Due to frustration effects, the ground-state structure is 
nontrivial for shifted parallel and perpendicular stripes. 
We now determine these ground states treating spins as 
classical. These ground states will be a necessary prereq- 
uisite for the subsequent spin- wave analysis. We continue 
to focus on the representative case p = 4. 

Depending on the values of the couplings A and p we 
find two different types of ground states. For a nearest 



neighbor exchange across a stripe in the range < A < A c 
(A c w 0.59 for shifted parallel and A c ~ 0.35 for per- 
pendicular stripes) the ground state has a canted planar 
topology up to a value ^ C (A) of the interlayer exchange 
(cf. Fig. 0). For p > p c {\) spins lock into a collinear 
texture. 
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7 


colli near 






canted 




planar 



FIG. 2: Classical ground-state phase diagrams for shifted par- 
allel and perpendicular stripes. For p < /ti c (A) the ground 
states show a canted planar spin pattern illustrated in Fig. Q 
For p > p c (X) the topology of the ground states changes into 
a collinear pattern where spins lying on top of each other 
are strictly antiparallel and nearest neighbors across a stripe 
are parallel. For A — > A c (Ac ~ 0.59 for shifted parallel and 
A c ~ 0.35 for perpendicular stripes) p c goes to infinity. Above 
Ac the ground states are always planar. 



To characterize these different phases, we start with 
the planar one. As already indicated above, the frustra- 
tion can lead to a canting of spins. The origin of the 
canting is easily understood. For p — the layers are 
decoupled and the sublattice magnetization in both lay- 
ers can have an arbitrary relative orientation. For small 
interlayer coupling p the spins start to cant starting from 
a configuration where spins lying on top of each other are 
perpendicular. Only in this case the interlayer couplings 
lead to an energy gain proportional to small canting an- 
gles while the intralayer couplings lead to an energy costs 
of second order in the canting angles. Such canted pla- 
nar ground states are illustrated in Fig. □ In Fig. the 
corresponding tilting angles are plotted for A = 0.1 as a 
function of p. The tilting angles increase monotonously 
in a way that spins lying on top of each other become 
increasingly antiparallel with increasing p. 

In the other phase, for p > /i c (A), the interlayer 
coupling p dominates the coupling A across the stripes 
and the topology of the ground state changes into a 
collinear configuration where the spins lying on top of 
each other are strictly antiparallel and nearest neighbor 
spins across a stripe are strictly parallel although they 
are antifcrromagnctically coupled. This configuration is 
stable against a canting of the spins because for small A 
the energy gain for A-bonds and the energy costs for p- 
bonds as well as the couplings within the domains would 
be quadratic in the tilting angles. Since this ground 
state has lost the antiphase-boundary character of the 
charge stripes it resembles a diluted antiferromagnet. 
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FIG. 3: Upper row: Energy per lattice site in units of JS as 
a function of fi for A = 0.1. For both stripe configurations the 
energies of the canted planar and the collinear spin pattern are 
plotted. The curves intersect at /x = fj, c where the topology 
of the ground states changes. Lower row: Relative values of 
the tilting angles of the spins in the planar configuration as a 
function of the interlayer coupling /j, for A = 0.1. 



This would lead to a static magnetic response at the an- 
tiferromagnetic wave vector in disagreement with exper- 
imental observations. Therefore, these collinear phases 
probably are not relevant for the magnetic properties of 
the cuprate compounds. 

For small values of A the phase boundary is approx- 
imately given by /x c (^) ~ 2 A for both stripe configura- 
tions (cf. Fig- O - In the limit A — > A c the critical value 
He goes to infinity. Above A c the ground states remain 
canted planar for all values of the interlayer coupling fi. 

Comparing the classical magnetic ground-state ener- 
gies for the two frustrated configurations, we find that - 
in contrast to the Coulomb energy - the exchange cou- 
pling favors perpendicular stripes over shifted parallel 
stripes. For this reason we retain perpendicular stripes 
in our consideration. 



III. SPIN- WAVE THEORY 

In this analytic part we derive general expressions for 
the magnon band structure and the spectral weight at 
zero temperature in a framework of linear spin-wave the- 
ory (for a review in the context of cuprates, see e.g. Ref. 
Ilq) . These expressions are evaluated numerically later on 
in Sec. lIVI for parallel, shifted parallel and perpendicular 
stripes and fixed stripe spacing p = 4. 



A. Holstein-Primakoff representation 

The ground-state analysis of the preceding section has 
made clear that spin waves now have to be introduced as 
excitation of a non-collinear ground state. However, our 
numerical calculation of the classical ground states have 
shown planar spin textures (here, a collinear texture is 
considered as a special subcase of a planar texture). 

In the following we consider a general planar ground 
state which can be captured by a vector field S Q (r) = 
{cos <^ a (r), sin 4> a (r), 0}, where the tilting angles of the 
spins obey the translational symmetry </> a (r) = <f) a (r+A) 
for an arbitrary magnetic lattice vector A = miAS 1 ' + 
7712 A^ 2 ) • For the spin textures displayed in Fig. ^ corre- 
sponding magnetic unit cells are given by A^ 1 ) = (4, 1) 
and A^ 2 ) — (0, 2) for parallel stripes and for shifted par- 
allel stripes, and by A^ = (8, 0) and = (0, 8) for 
perpendicular stripes. 

To study the quantum fluctuation around the classical 
ground state we rotate all spins by their planar angles 
(j> a (r) according to 

SI (r) = S*(r)cos<k(r)-S»(r) S in0 Q (r), (4a) 
S*(r) = S*(r)sin0 Q (r) + ^(r)cos^(r), (4b) 
Sl(v) = Sl(v), (4c) 

such that S(r) has a classical ferromagnetic ground state 
S(r) = S'jFO, 0}. In the transformed spin basis we in- 
troduce Holstein-Primakoff (HP) bosons in the standard 
way (using S± = S v ±iS z ) 7 



(5a) 
(5b) 
(5c) 



S+(r) = ^j2S-h r , a b r , a , 

S*(r) = -n T , a +S, 
and obtain the spin-wave Hamiltonian 

r,r' a, a' 

(r,r') + &L&L']}>( 6 ) 

where the functions / and g are defined by 

/ Q , Q <(r,r') = i [J a (r,r')6 a , a ' + ^JS r , r ' (1 - 5 a , a >)] 
x [A Q , Q ,(r,r') + 1] 

r f 

-yJSr,r'S a ,a> ~ S a,a> ) A a ,a> (r, r') 

a' 

(7a) 

g a . a >(r,r') = - [J a (r, t')5 a ^ + /J,J5 rtT > (1 - £ a ,a')] 

x [A Q , Q -(r,r')-l] (7b) 
A a>a /(r,r') = cos[0 a (r) - a /(r')] . (7c) 



5 



To diagonalize the Hamiltonian, we Fourier transform the 
bosonic operators via b a (r) — f, exp(ikr)6 ct (k), where 
J, = (27r) -2 J cPk and the k integrals run over the Bril- 
louin zone of the square lattice with an area (2-7r/a) 2 . 
Following our calculations for the monolayer system-^ 
we decompose a square lattice vector r into a magnetic 
lattice vector A and a decoration vector a (r = A + a) . 
The number of vectors a is denoted by n (the area of the 
magnetic unit cell). In momentum space, the reciprocal 
magnetic basis QW , i = 1,2, spans the corresponding 
magnetic Brillouin zone (BZ). Wave vectors k can be 
uniquely decomposed into k = Q + q with q G BZ and 
Q = toiQW +m 2 Q (2) . Within the Brillouin zone of the 
square lattice there are n vectors Q which we denote by 
Q„. Using these decompositions we rewrite the spin- wave 
Hamiltonian as 



n 



J EE F ^, l /'a'(q)[^,q + Q„ V,q+Q„, 
v.v' a.a' 



-6a,-q-Q„&L 



q-Q„ 



+\ J E E G -x^q)[^, q+Q ,X',-q-Q,v 

*^ a, a' 

+&a,-q-Qj,&a',q+Q„/]j (8) 



where 

Fva^v' a' (q) 



s 



-EE/«.«'( a + A ' a ') 



n 

A a.a' 

x cos [qA + q(a - a') + Q^a - Q,ya(; 



is essentially the Fourier transform of /, 
g 

-fa,c*'(Qv + <l,Qu> +q') =S(q + q')F ua y^(q L ). (10) 
n 

Analogous expressions relate G to g. The Hamiltonian 
||SJ has exactly the same structure as in the monolayer 
case [compare Eq. 8 in Ref. and can be diagonal- 
ized by a Bogoliubov transformation in an analogous way. 
The final diagonal form is given by 



n = f^J^(q)( [ bt ( ( q )b J (< l ) + lY 



(11) 



where the squared energies w 2 are eigenvalues of the her- 
mitian matrix M- 1 / 2 KM~ 1 / 2 . Thereby 1VU 1 = F — G 
denotes the inverse mass matrix and K = F + G the 
coupling matrix. 



B. Structure Factor 

We now proceed to calculate the inelastic zero- 
temperature structure factor for even and odd excitations 



S£(k,u;) 



E E 

F j=x,y,z 



|<F|^(k)±^(k)|0)| 2 
xS(cj-lj f ). (12) 



Here, |0) denotes the ground state (magnon vacuum) 
characterized by 6 7 (q)|0) = and we consider only 
single- magnon final states |F) = 6Uq)|0) with excita- 
tion energy uj-p := Ep — Eq. k = (k x ,k y ) denotes 
the in-plane wave- vector, odd excitations correspond to 
k~ = (2n + l)n/d [L~ = (2n + l)c/(2d) in reciprocal 
lattice units], even ones to fc^ = 2m: /d (L + = nc/d), 
where d is the distance of the two layers within the or- 
thorhombic unit cell. For YBCO with d w 3.34 A and 
c w 11.7 A the corresponding values for even and odd 
modes arc L~ rj 1.75, 5.25 and L + » 0, 3.5. 

Expressing the spin operators by the final bosonic 
operators 6 7 (q) it is straightforward to calculate the 
structure factor. Using a pseudo-Dirac notation and 
denoting the 2n-dimensional cartesian basis by \v, a)) 
(y = l,...N,a — 1,2) and the orthonormal eigenba- 
sis of M _1 / 2 KM -1 / 2 by I7)), the structure factor can be 
rewritten in a compact form, 

S^q + Q^w) = 5^5 7 ± (q + Q,) ( 5(a;-a; 7 (q)),(13a) 

7 

^(q + Q.) = \ i^±\™- 1/2 \l)) 
X=C,S,w 7 1 K 

xi-«7|M- 1 /»X|i/,±>, (13b) 
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where we have defined \u, ±)) = (l/y2)[\u, 1)) ± \u, 
and introduced the matrices S and C according to 



^vol^v' a' 



^vol^v' a' 



-^Vsin0 a (a)e i ^-^') a , (14a) 

a 

1 

-8 aa > Vcos0 Q (a)e l(Q -- Q -' )a . (14b) 

n — ' 



IV. RESULTS 

We now evaluate the the magnon dispersion and the 
inelastic structure factor for even and odd excitations 
numerically. From a comparison of our findings for the 
monolayer system to neutron scattering data for the 
cuprate compounds we founcU^ the coupling \J across 
a stripe to be about one order of magnitude smaller than 
the nearest neighbor coupling J within the domains. For 
the coupling \iJ between the layers a value /1 » 0.08 is 
reported^ for antiferromagnetic YBCO in the absence of 
stripes. Therefore in the stripe system the couplings A 
and n can be assumed to be of the same order. In the 
following we keep the value of A fixed and discuss the ef- 
fects of increasing fi starting from the case of decoupled 
layers (/i = 0) where the band structure of the mono- 
layer systemic should be recovered. In this parameter 
regime the classical ground states for shifted parallel and 
perpendicular charge stripes show the canted planar tex- 
ture and the antiphase domain boundary character of the 
charge stripe is weakened by the interlayer coupling but 
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still pronounced. Finally we shortly present the excita- 
tion spectra for shifted parallel and perpendicular stripes 
for parameters belonging to the collinear ground state 
regime. 

In the case of decoupled layers (fi = 0) the results of 
the monolayer system are trivially recovered. Since the 
two layers are uncorrelated, the structure factor does not 
depend on the L component of the wave vector. For 
parallel stripes (with or without a relative shift of the 
stripes) where the charge modulation is unidirectional 
with Qf 1 = Qf 1 = (1/4,0) we just obtain an additional 
twofold degeneracy of each of the three bands due to the 
equivalence of the two layers. Therefore the degenera- 
tion of the bands is fourfold since in the monolayer case 
each band is twofold degenerated due to the equivalence 
of the two sublatticesii^ The lowest, acoustical band has 
zeros at the magnetic superstructure vectors which are 
located at (j/4,0) and (j/4 + 1/8,1/2), j = 0, ...,3, 
within the Brillouin zone of the square lattice (we choose 
< H, K < 1). The spectral weight is concentrated 
near the lowest harmonic incommensurate wave vectors 
Q = (l/2± 1/8, 1/2). With increasing energy the incom- 
mensurability decreases and the branches of the acoustic 
magnon band close at the antiferromagnetic wave vector 
(1/2, 1/2) and an energy which we associate with the 
7r-resonance. Along the (H,l/2) direction the acoustic 
band is gapped to the overlying optical magnon band (see 
upper left panels in Figs. and I^Jl. Along the orthogonal 
direction (1/2, K), one optical band has vanishing spec- 
tral weight and only two bands are visible (see middle-left 
panels in Figs. 0] and [5J. 

In twinned samples with stripe domains oriented or- 
thogonal to each other, a scan along the (H, 1/2) di- 
rection results in the superposition of the signals ob- 
tained from scans in directions (H, 1/2) and (1/2, H) of 
a single-domain sample. For domains of equal size, one 
thus obtains an apparent symmetry (H, K) <-> (K, H) 
and a fourfold pattern of the static incommensurate 
wave vectors located at Q = (1/2 ± 1/8, 1/2) and Q = 
(1/2, 1/2 ± 1/8) also for (shifted) parallel stripes. In 
Figs.0]and[5l the panels in the third row are just obtained 
by superimposing the panels of the first and second row. 
Since the acoustic band of the monolayer system has a 
saddle point at the antiferromagnetic wave vector, the re- 
sulting band structure is x-shaped in the vicinity of the 
7r-resonance energy. 

The configuration of hole stripes lying perpendicular to 
each other corresponds to charge modulation wave vec- 
tors Qf = (1/4,0) and Qf = (0,1/4). For decoupled 
layers, the resulting band structure contains the bands 
of the monolayer system and the same bands rotated by 
90 degrees leading to the symmetry ui(H, K) — uj(K, H) 
and therefore to a fourfold pattern of the static incom- 
mensurate wave vectors located at Q = (1/2 ± 1/8, 1/2) 
and Q = (1/2, 1/2 ± 1/8). Thus, for /1 = 0, the structure 
factor is identical for perpendicular stripes and twinned 
parallel stripes (left lower panel in Figs. 01 and EJ. 

With increasing interlayer coupling /i the bands start 
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FIG. 4: Band structure and spectral weight along the 
(H, 0.5, L*) and (0.5, H, L*) directions for parallel stripes ly- 
ing on top of each other and couplings A = 0.15 and /1 = 0, 
0.08. The last row shows the band structure of a twinned 
sample (see text). L + corresponds to even, L~ to odd excita- 
tions. Darker and larger points correspond to a larger weight 
of the inelastic structure factor. 



to split with different distributions of the spectral weights 
in the odd and even channel (cf. Figs. 01 03 and|S}. For 
parallel and shifted parallel stripes the Hamiltonian is 
invariant under the reflection z — + — z combined with a 
translation. This implies that the magnon states - mod- 
ulo a phase factor which does not enter the structure 
factor - have a well defined parity with respect to an 
exchange of both layers. As a consequence, nondegener- 
ate bands are visible only either in the even or the odd 
channel. 

Nevertheless the excitation spectra of the two parallel 
stripe configurations deviate significantly, e.g. the even 
excitations are gapped for parallel stripes whereas for 
shifted parallel stripes the intensity of even excitations 
is only reduced at low energies (cf. middle columns of 
Figs. 01 [SJ. For stripes on top of each other each band - 
which is fourfold degenerate at \x — - splits up into 
twofold degenerate bands which have identical parity. 
For shifted stripes each band splits up into three bands. 
One of them is twofold degenerate and both subbands 
are of opposite parity. Therefore this degenerated band 
is visible in both channels (cf. Fig. [SJ . 

For perpendicular stripes the symmetry is more com- 
plicated. The Hamiltonian is invariant under a reflec- 
tion z — > — z in combination with a 90°-rotation along 
the z axis. Since this rotation mixes different wave vec- 
tors, almost all eigenstates do not have a well defined 
parity and will be partially visible in the odd and even 
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FIG. 5: Band structure and spectral weight along the 
(H, 0.5, L ) and (0.5, H,L ) directions for shifted parallel 
stripes and couplings A = 0.15 and /i = 0,0.08. The last 
row shows the resulting band structure of a twinned sample. 




FIG. 6: Band structure and spectral weight for even (left 
panel) and odd (right panel) excitations along (H, 0.5, L*) 
direction for perpendicular stripes with couplings A = 0.15 
across the stripes and interlayer couplings /j, = 0.04 (upper 
row) and fi = 0.08 (lower row). 



channel. The exception are modes at particular wave 
vectors such as the antiferromagnetic wave vector which 
are mapped onto themselves (modulo a reciprocal lattice 
vector). Only there the excitations can be classified due 
to their symmetry. Like for the shifted parallel stripes 
the excitations are not gapped in the even channel (cf. 
Fig.EJ. 

We now focus on the band splitting and the distribu- 
tion of the spectral weights of even and odd excitations 
at the antiferromagnetic wave-vector (1/2,1/2). With 
increasing interlayer coupling fi, the resonance energy 
ui n splits up into two different energies uj~ and w+ for 
centered-parallel stripes and into three energies w~, 
and w+ for the other stripe configurations as schemat- 
ically illustrated in Fig. [7| It is common to all stripe 
configurations that uj~ has a finite spectral weight only 
in the odd channel, whereas uj^ has a finite weight only 
in the even channel. For shifted parallel and perpendic- 
ular stripes, in both channels a finite intensity is found 
at the intermediate energy u®. This intensity is however 
smaller than at u>^. 

The splitting of the resonance energy for shifted par- 
allel and perpendicular stripes looks quite similar. lu~ 
and a>+ are almost equidistant to the intermediate en- 
ergy w° which increases only slightly with fj, (cf. Fig.|SJ). 
For small couplings the splitting is quadratic in /i. For 
centered-parallel stripes the splitting looks different, oj+ 
increases almost linearly with fi whereas oj~ is almost 
independent of the interlayer coupling. 

Finally, we calculate the band structures for shifted 
parallel and perpendicular stripes for couplings A < A c 
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FIG. 7: Schematic illustration of the band splitting in the 
vicinity of the antiferromagnetic wave vector (1/2, 1/2) along 
the (H, 1/2) direction. In the cases of parallel stripes the 
band structures for twinned samples are shown. Even and 
odd bands are gathered together. 



and n > /U C (A) where the ground-states are collinear and 
the charge stripes lose their anti-phase domain boundary 
character. We implicitly assume that fj, is not too large, 
otherwise spins on top of each other dimerize and lose 
their magnetization. In this regime the magnetic fluctu- 
ations are drastically changed. For both stripe orienta- 
tions, the odd channel now has a static signal at the an- 
tiferromagnetic wave vector, whereas in the even channel 
the spectral weight is concentrated at incommensurate 
positions (1/2 ± 1/4, 1/2) (cf. Fig.EJ. For perpendicu- 
lar stripes we also find small intensity at this positions 
in the odd channel. The incommensurability is doubled 
compared to the regime of canted planar ground-states 
reflecting that the charge stripes do not act like anti- 
phase domain boundaries in the regime of strongly cou- 
pled layers. In the even channel the intensity at the an- 
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FIG. 8: Splitting of the resonance energy as a function of 
the interlayer coupling /1 for A = 0.15. In the odd channel 
(L = L~) the spectral weight is concentrated at oj^ and no 
intensity is found at , in the even channel (L = L + ) no 
excitations at ui~ are observable and the spectral weight is 
concentrated at ui%. For shifted and perpendicular stripes in 
both channels a small intensity is found at the intermediate 
energy u°. 



tiferromagnetic wave-vector is peaked at an energy ui n 
which increases with the interlayer coupling /1 and is ap- 
proximately the same for both stripe configurations. 
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FIG. 9: Band structure in the collinear regime fi > /i c (A) 
along the direction (H, 1/2, L ) for shifted parallel and per- 
pendicular stripes with A = 0.07 and fj, = 0.50. 



V. DISCUSSION 

In this section we compare our results to neu- 
tron scattering data for the bilayer high-T c compound 
YBa 2 Cu 3 06+ x - We wish to stress that - because of the 
simplifications assumed in our model - it is not our goal 
to obtain a quantitative agreement. Rather we wish to 
draw a qualitative comparison in order to fortify the hy- 



pothesis that the stripe picture is a suitable approach to 
describe spin fluctuations. Furthermore, we hope that a 
comparison of future experimental data with our calcu- 
lations will help to identify the realized stripe configura- 
tion. 

Since a spin gap with an energy w ga p - e.g. due 
to Cooper-pair formation - is not incorporated in our 
model, the results apply only to energies above w gap 
where the magnon dispersion is not masked by the super- 
conducting condensate. In particular in the underdoped 
regime where w gap decreases with the doping level, the 
calculated spectral features become visible over an in- 
creasing energy range. Our calculations are restricted 
to zero temperature. Therefore, a comparison can also 
be made only to experiments performed at temperatures 
well below the superconducting transition temperature. 

ExperimentsiZii£ii£ in (partially) detwinned YBCO 
provide evidence for unidirectional order, i.e., that a four- 
fold pattern of incommensurate peaks near the antifer- 
romagnetic wave vector Icaf = (1/2,1/2) results only 
from the twinning. The stripes seem to be parallel and 
oriented along the direction of the oxygen chains in the 
adjacent planes. This immediately speaks against the 
scenario of perpendicular stripes for which detwinning 
would not affect the fourfold symmetry. 

We briefly recall some neutron scattering measure- 
ments on YBa2Cu306+ x which provide insight into the 
incommensurability and the 7r resonance over a wide dop- 
ing and temperature range. It was controversial for quite 
some time whether both phenomena would exist above 
T c until in underdoped materials the incommensurability 
was found also above T c m& Likewise, the appearance of 
the magnetic resonance was found above T c , occurring 
together with the pseudogap at a temperature T* > T c 
determined from transport and nuclear resonance^i Al- 
though the 7r-resonance persists as a well defined feature 
also in the normal state above T c , its intensity can be 
reduced significantly at T c i& For near optimally doped 
compounds, the resonance is not detectable in the normal 
phased since T* almost coincides with T c . Dai el alm& 
concluded that the resonance exists above T c for x < 0.8 
and that incommensurate spin fluctuations appear in the 
normal state for x < 0.6. Arai et alw^ also observed 
incommensurate fluctuations in the normal state for a 
sample with an oxygen concentration of x = 0.7. Thus, 
superconductivity is not a prerequisite for incommensu- 
rability and 7T resonance in bilayer compounds as well as 
in monolayer compounds. 

For underdoped YBCO with various oxygen concentra- 
tions, the experimentally observed spin dynamics data 
(see Table HJ look qualitatively very similar. There is 
a systematic increase of the incommensurability and of 
the 7T resonance frequency with doping, which is consis- 
tent with our model. We have shown this recently for a 
monolayer model^ The bilayer stripe model shares this 
feature and therefore we focus in this paper exclusively 
on specific bilayer features. 

Experimentally, constant energy scans slightly above 
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TABLE I: Spin dynamics data for YBa2Cu30s+x for various 
oxygen concentrations x characterized by the critical temper- 
ature T c , incommensurability 8, corresponding stripe period 
p, the resonance energy uj~ observed in the odd channel, and 



the gap in the odd channel along (H, 1/2, L~) show 
a broad intensity peak at kAF, before incommensurate 
scattering sets in and the data can be compared to our 
model. The intensity shows magnetic peaks at a dis- 
tance 5k(u>) away from k^F- The incommensurability 6 
is determined by extrapolating 5k(uj) to lo = and it is 
connected to the stripe spacing p through 6 — l/(2p). 
The incommensurate peaks are best defined if the stripe 
spacing is nearly a multiple of the lattice spacing (integer 
p) since the stripes are stabilized by the latticed 

The three stripe configurations examined for our model 
are not equivalent in their low-energy behavior. For (un- 
shifted) parallel stripes (see Fig.0J, an incommensurabil- 
ity is visible at low energies only in the odd channel since 
the even channel has a relatively large gap not related to 
superconductivity. In contrast, for shifted parallel and 
perpendicular stripes the even channel shows incommen- 
surate response down to the superconducting gap. Ex- 
perimental evidence2ii2Si2& for a large gap in the even 
channel (well above the resonance energy in odd channel) 
therefore favors the configuration with unshifted parallel 
stripes. 

With increasing energy, the separation 5k(tu) of the in- 
commensurate peaks decreases and the branches close at 
kAF at certain energies w,. Depending on the stripe con- 
figuration, there are two or three such energies, compare 
Fig. [7| According to our model, an energy scan of the 
odd channel at kAF would show a first resonance at the 
intersection with the lowest magnon band at lu~ which 
we identify with the resonance frequency&21 For shifted 
parallel and perpendicular stripes, a second line at oj® 
contributes to the odd channel. It has significantly less 
weight and is separated from the first one by only a small 
energy splitting (of the order of a few meV) which would 
be hard to be resolved experimentally. 

In a similar way, the even channel has a resonance at 
an energy u>+ > uj~ , and for shifted parallel an perpen- 
dicular stripes also a weaker resonance at an intermedi- 



ate frequency wJJ (cf. Fig. EJ) . Experimentallyfi2i2L2& a 
strong oscillatory dependence of the scattering intensity 
on L shows that the resonance frequencies in the odd and 
even channel are well separated. Energy scans at the an- 
tiferromagnetic wave vector show peaks at tu~ in the odd 
channel and oj+ in the even channel, no peak at the inter- 
mediate energy which should be visible in both chan- 
nels is resolved^ This again favors unshifted parallel 
stripes, which (in contrast to shifted parallel and perpen- 
dicular stripes) have no shared resonance frequency 
Although we restricted our comparison to experiments 
on underdoped samples, overdoped compounds also show 
two distinct resonance modes of opposite symmetry, 29 
which could be identified with uj~ and u>^ . 

From a comparison of the band splitting Alu^ = lu+ — 
uj~ to experimental values (cf. Tab. we can estimate 
the strength of the interlayer coupling /i. For A = 0.15 
we find fi w 0.02 — 0.06 almost independent of the stripe 
configuration. This value is reasonable since the effective 
coupling fi in the stripe system should be slightly reduced 
compared to the undoped case where a value of \x « 0.08 
is reported^ 

Above lOtt the response is found to become incommen- 
surate again with increasing separation Sk(uj). The mo- 
mentum width is larger and the intensity is weaker than 
below LQ-K- Overall, the dispersion is "x-shaped". As 
pointed out in Sec. IIVI such a shape appears basically 
for every non-unidirectional stripe configuration, for par- 
allel stripes in twinned crystals as well as in perpendic- 
ular stripes. The x-shape has been observed explicitly 
in Refs. l7 U24ll25ll26l It would be interesting to verify in 
detwinned samples that the relative intensities of the up- 
per and lower branches of the x-shape are related to the 
population ratio of the twin domains. 

In conclusion, we have calculated the bilayer effects in 
the magnetic excitation spectrum in striped states. As a 
generic feature of the stripe model we find x-shaped dis- 
persion relation in the vicinity of the it resonance, which 
is consistent with experimental data. We have obtained 
a bilayer splitting of single-layer bands into two or three 
bilayer bands. From the three stripe configurations stud- 
ied, the unshifted parallel case overall is most consistent 
with neutron scattering data, although it seems to be 
energetically unfavorable at first sight. 
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